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Abstract 
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The masses of the neutral CP-even Higgs bosons in the Minimal Super- 
symmetric Standard Model (MSSM) are predicted on the basis of explicit 
Feynman-diagrammatic calculations. The results, containing the complete 
diagrammatic one-loop corrections, the leading two- loop corrections of (D(aa s 
and further improvements taking into account leading electroweak two-loop 
and higher-order QCD contributions, are discussed and compared with re- 
sults obtained by renormalization group calculations. Good agreement is 
found in the case of vanishing mixing in the scalar top sector, while sizable 
deviations occur if scalar top mixing is taken into account. By means of 
a Taylor expansion a compact approximation formula for the mass of the 
lightest Higgs boson, m^, is derived. The quality of the approximation in 
comparison with the full result is analyzed. 
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The masses of the neutral CP-even Higgs bosons in the Minimal Su- 
persymmetric Standard Model (MSSM) are predicted on the basis of ex- 
plicit Feynman-diagrammatic calculations. The results, containing the 
complete diagrammatic one-loop corrections, the leading two-loop correc- 
tions of 0(aa s ) and further improvements taking into account leading 
electroweak two-loop and higher-order QCD contributions, are discussed 
and compared with results obtained by renormalization group calculations. 
Good agreement is found in the case of vanishing mixing in the scalar top 
sector, while sizable deviations occur if scalar top mixing is taken into ac- 
count. By means of a Taylor expansion a compact approximation formula 
for the mass of the lightest Higgs boson, mh, is derived. The quality of the 
approximation in comparison with the full result is analyzed. 

PACS numbers: 11.30.Pb, 12.38.Bx, 14.80. Cp 

1. Introduction 

The search for the lightest Higgs boson provides a direct and very strin- 
gent test of Supersymmetry (SUSY) and is one of the main goals at the 
present and the next generation of colliders. A precise prediction of its 
mass, mh, is inevitable for determining the discovery and exclusion poten- 
tial of LEP2 and the upgraded Tevatron in this search and for analyzing the 
accessible MSSM parameter space. If the MSSM Higgs boson exists, it will 
be detectable at the LHC and a future linear collider (LC), and its mass 
will be measured at these machines with high precision. The comparison 
of the MSSM prediction with the experimental value of m-h will then allow 

* Talk given by G. Weiglein at the Cracow Epiphany Conference on Electron-Positron 
Colliders, Cracow, January 5-10, 1999. 
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a very sensitive test of the model. A precise knowledge of the mass of the 
heavier CP-even Higgs boson, mn, will be important for resolving the mass 
splitting between the CP-even and -odd Higgs-boson masses. 

The mass of the lightest Higgs boson in the MSSM is restricted at the 
tree level to be smaller than the Z-boson mass, Mz- The dominant one-loop 
corrections arise from the top and scalar-top sector via terms of the form 
Gf^mf ln(m{ nif /mj ) [1]. These results have been improved by perform- 
ing a complete one- loop calculation in the on-shell scheme [2-4]. Beyond 
one-loop order renormalization group (RG) methods have been applied in 
order to obtain leading logarithmic higher-order contributions [5-9]. Fur- 
thermore the leading two-loop QCD corrections have been calculated in the 
effective potential method [10,11]. Phenomenological analyses for the neu- 
tral CP-even Higgs-boson masses have until recently been based either on 
RG improved one- loop calculations [6, 7, 9] or on the complete Feynman- 
diagrammatic one-loop on-shell result [2-4] . The numerical results of these 
approaches however differ by up to 20 GeV in m^. 

Recently the Feynman-diagrammatic result for the dominant two-loop 
contributions of 0(aa s ) to the masses of the neutral CP-even Higgs bosons 
has become available [12]. By combining these contributions with the com- 
plete one- loop on-shell result [3], the currently most precise result for mj, 
based on diagrammatic calculations is obtained [13,14]. It has been imple- 
mented into a Fortran program called FeynHiggs [15]. In the present paper 
the new Feynman-diagrammatic results are briefly summarized and com- 
pared with the results obtained by RG methods. Furthermore a compact 
analytical approximation formula [16] is discussed, which is derived from 
the full diagrammatic result by means of a Taylor expansion. 



2. Diagrammatic two-loop calculation of the masses of the 
neutral CP-even Higgs bosons 

The MSSM Higgs sector can be described with the help of two para- 
meters: tan/3 = V2/V1, the ratio of the two vacuum expectation values, and 
Ma, the mass of the CP-odd Higgs boson. The tree-level predictions for 
the masses and ttt-h of the neutral CP-even Higgs bosons h and H are 
determined by diagonalizing the tree-level mass matrix given in terms of the 
current eigenstates <pi and fa- ln the Feynman-diagrammatic approach the 
higher-order corrected Higgs-boson masses are derived by determining the 
poles of the h, i^-propagator matrix whose inverse is given by 
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where the S denote the renormalized Higgs-boson self-energies, which can 
be decomposed according to 

t s = SW + + . . . , s = h,H,hH, (2) 

into the contributions at one-loop order, two-loop order etc. 

For the one-loop contributions to these self-energies, si 1 ^ 2 ), we take 
the result of the complete one- loop on-shell calculation of Ref. [3]. The 
agreement with the result obtained in Ref. [2] is better than 1 GeV for 
almost the whole MSSM parameter space. 

The leading two-loop corrections, Tr s (0), have been obtained in Refs. 
[12-14] by calculating the 0(aa s ) contribution of the t,t-sector to the 
renormalized Higgs-boson self-energies at zero external momentum from the 
Yukawa part of the theory. The calculation has been performed in the on- 
shell scheme. It involves a two-loop renormalization in the Higgs sector and 
a one-loop renormalization in the scalar top sector of the MSSM. The cal- 
culations have been performed using Dimensional Reduction (DRED) [17], 
which is necessary in order to preserve the relevant SUSY relations. In 
deriving these results, use has been made of the computer-algebra pro- 
grams FeynArts [18] (in which the relevant part of the MSSM has been 
implemented) for generating the Feynman amplitudes, and TwoCalc [19] 
for evaluating the two-loop diagrams and counterterm contributions. 

The results for the corrections in 0{aa s ) are given in terms of the SUSY 
parameters tan/3, Ma, fJ,, nig, m^, m^ 2 , and 9^, where [i denotes the Higgs- 
mixing parameter and nig the mass of the gluino. The mass eigenstates t\, t<i 
and the mixing angle 9~ t in the scalar top sector are derived by diagonalizing 
the mass matrix of the scalar top quarks given in the basis of the current 
eigenstates it, £r. The non-diagonal entry in the scalar quark mass matrix 
is proportional to the mass of the quark and reads for the t-mass matrix 
mtM[ R = mt(A t — /Ltcot/3), where we have adopted the conventions as 
in Ref. [14]. Due to the large value of m t these mixing effects are in general 
non-negligible. 

Inserting the contributions in 0(a) and 0(aa s ) into eq. (1) and deter- 
mining the poles of the h, ff-propagator matrix yields the prediction for the 
masses of the neutral CP-even Higgs bosons. We have implemented two 
further corrections beyond (D(aa s ) into the prediction for m^: The first 
correction concerns leading QCD corrections beyond two-loop order, taken 
into account by using the MS top-quark mass 

m t = m t (m t ) ^m t / (l + ^-a s (m t )\ (3) 



for the two- loop contributions instead of the pole mass m t . The second one 
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is the leading two-loop Yukawa correction of 0{G 2 l m^), taken over from the 
result obtained by RG methods [7, 20] . 

The results described above have been implemented into the Fortran 
program FeynHiggs [15], which needs about 0.5 seconds for the evaluation 
of rah, ran on a Sigma station (Alpha CPU, 600 MHz) for one set of pa- 
rameters. As an additional constraint (besides the experimental bounds) 
on the squark masses, the program also evaluates the contribution to Ap 
arising from i/6-loops up to 0(aa s ) [21]. A value of Ap outside the preferred 
region of Ap SUSY < 1 • 10~ 3 [22] indicates experimentally disfavored t- and 
6-masses. The program FeynHiggs is available via the WWW page 
http : //www-itp .physik.uni-karlsruhe .de/f eynhiggs . 

3. Numerical Results 

For the numerical evaluation we have chosen two values for tan (3, which 
are favored by GUT scenarios [23]: tan/3 = 1.6 for the SU(5) scenario 
and tan/3 = 40 for the SO(10) scenario. Other parameters are M% = 
91.187 GeV,M w = 80.39 GeV, = 1.16639 10" 5 GeV- 2 ,a s (m t ) = 0.1095, 
and m t = 175 GeV, if not otherwise indicated. Further parameters are 
M A , m g , p, and the SU(2) soft SUSY-breaking parameter M{= M 2 ). The 
other gaugino mass parameter, Mi, is fixed via the GUT relation M\ = 
(5 fi w)/(3cw) M. In the figures below we have chosen nig = M$ L = Mf R for 
the diagonal entries in the scalar top mass matrix. 
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Fig. 1. One- and two- loop results for the mass of the lightest Higgs boson mj, as a 
function of M^ R /rriq for two values of tan/3. 

Fig. 1 shows the result for rrih obtained from the diagrammatic two-loop 
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calculation as a function of M f /rriq, where rriq is fixed to 1000 GeV. The 
two-loop contributions give rise to a large reduction of the one-loop on-shell 
result by up to 20 GeV. A minimum in the prediction for occurs around 
M^ R /rriq = 0, which corresponds to the case of no mixing in the t-sector. A 
maximum in the two-loop result for mj, is reached for about |M t /rriq\ ~ 2, 
this case we refer to as 'maximal mixing'. In the two-loop result the maxima 
are shifted compared to their one-loop values of about \M^ R /rriq\ ~ 2.4. 
Varying tan/? around the value tan/3 = 1.6 leads to a relatively large effect 
in rrih, while the effect of varying tan/3 around tan/3 = 40 is marginal. 
Different values of the gluino mass, nig, in the two-loop contribution affect 
the prediction for by up to ±2 GeV in the maximal-mixing scenario, 
while the effect is negligible in the no-mixing scenario. Varying M, which 
enters via the non- leading one-loop contributions, changes the value of 
by ±1.5 GeV. A more detailed analysis of the dependence of our results on 
the different SUSY parameters has been performed in Ref. [14] . 
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Fig. 2. The two-loop result for mi, as a function of m q in the no-mixing and the 
maximal-mixing case. The point marked by an arrow indicates the prospect for 
the experimental precision reached by a future linear collider in the determination 
of rrih and m q for the hypothetical values mj, = 115 GeV and m q = 400 GeV. 

If the lightest Higgs boson and Supersymmetric particles will be found 
at the next generation of colliders, the experimental value of will be 
measured with high accuracy and also the possible range of the SUSY scale 
rriq will in this case be constrained to a small interval. At a high-luminosity 
LC the prospect for the accuracy obtainable for these parameters is Am^ = 
0.05 GeV and Arriq = 0.1%. In Fig. 2 the two-loop result for is shown 
as a function of rriq in the no-mixing and the maximal-mixing case. The 
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parameter space in the (m^, nig) plane corresponding to the accuracy in mj, 
and rriq at the LC is indicated in the plot for the hypothetical central values 
m-h = 115 GeV and m 9 - = 400 GeV. As can be seen from the plot, a precision 
determination of and nig will provide a very sensitive consistency test 
of the model. 
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Fig. 3. The maximally possible value for m,h as a function of tan/3 for m q — 
1000 GeV and three different values of m t . 



In order to determine the maximally possible value for mj, within the 
MSSM as a function of tan f3, we have performed a parameter scan in which 
rrig, M, fj,, Ma and M^ R have been varied for three values of mt and fixed 
values of nig and tan (3. Fig. 3 shows the maximal Higgs-boson mass value 
in the range tan (3 < 5 for rag = 1000 GeV (in Fig. 3 the choice M = \x = 
has been made for simplicity; the change in mj, when these parameters 
are chosen at their experimental lower bounds is negligible). The upper 
bound is shown for the current experimental value of the top-quark mass, 
mt = 173.8 GeV, and for values which are higher by one and two stan- 
dard deviations, respectively. Our results confirm that for the scenario with 
tan/3 = 1.6 practically the whole parameter space of the MSSM can be 
covered at LEP2. For slightly larger tan/3 and maximal mixing, however, 
some parameter space remains in which the Higgs boson could escape the 
detection at LEP2. For tan f3 = 40, on the other hand, the prediction for 

is at the edge of the LEP2 range even in the no-mixing case. The full 
exploration of the MSSM parameter space for the scenario with large tan/3 
will be a challenge for the upgraded Tevatron, the LHC, and the LC. 
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4. Numerical comparison with the RG approach 

We now turn to the comparison of our diagrammatic results with the pre- 
dictions obtained via RG methods. The upper plot of Fig. 4 shows the pre- 
diction for mj, as a function of M^ R /rriq, corresponding to our diagrammatic 
result and to the result obtained by RG methods [8] . In the no-mixing case 
the diagrammatic result agrees well with the RG result. For non- vanishing 
t-mixing sizable deviations between the diagrammatic and the RG results 
occur, which can reach 5 GeV for moderate mixing and become very large 
for large values of \M^ R /rriq\. As already stressed above, the maximal value 
for mj, in the diagrammatic approach is reached for \M^ R /m,q\ « 2, whereas 
the RG results have a maximum at \M^ R /rriq\ pa 2.4, i.e. at the one-loop 
value. Varying the value of mg in our result shifts the diagrammatic re- 
sult relative to the RG result (which does not contain the gluino mass as a 
parameter) within ±2 GeV in the region of large mixing. 

In the upper plot of Fig. 4 the results of our diagrammatic on-shell 
calculation and the RG methods have been compared in terms of the pa- 
rameters Mf , and M^ R of the i-mixing matrix. However, since the 
two approaches rely on different renormalization schemes, the meaning of 
these (non-observable) parameters is not precisely the same in the two ap- 
proaches starting from two-loop order. In order to compare results obtained 
by different approaches making use of different renormalization schemes, we 
find it preferable to compare predictions for physical observables in terms 
of other observables (instead of unphysical parameters). As a step into this 
direction we compare in the lower plot of Fig. 4 the diagrammatic results 
and the RG results as a function of the physical mass m~ t2 and with the mass 
difference Am ( - = m^ 2 — mi and the mixing angle 6~ t as parameters. In the 
context of the RG approach the running t-masses, derived from the i-mass 
matrix, are considered as an approximation for the physical masses. As in 
the comparison performed above in terms of unphysical parameters, in the 
lower plot of Fig. 4 very good agreement is found between the results of the 
two approaches in the case of vanishing t-mixing. For the maximal mixing 
angle 6\ = — tt/4 (and Am ( - = 340 GeV, for which the maximal Higgs-boson 
masses are achieved), however, the diagrammatic result yields values for mj, 
which are higher by about 5 GeV. 

The upper bound on m\ for a certain value of tan [3 derived from our 
diagrammatic results is thus higher in the low tan (3 region by about 5 GeV 
than the upper bound derived previously from the RG results. As a result, 
we find that the tan/3 region which can fully be covered at LEP2 and the 
upgraded Tevatron is significantly reduced compared to previous studies. 



8 



140 



120 



100 

> 

CD 

C5 80 



60 



40 




20 
-3.0 





110 




100 




90 




80 


> 

CD 


70 


U 






60 


E 






50 




40 




30 




20 



100 



■■ 1 000 GeV, m - = 500 GeV, n = -m-, M = m -, M n = 1 000 GeV 



2-loop RG 

2-loop diag., incl. all ref. 



-2.0 



-1.0 



M 



0.0 

'/rru 



1.0 



2.0 



3.0 




= 1 000 GeV, m~ = 500 GeV, tanp = 1 .6 
M 2 = 400 GeV, n = -200 GeV 



e f = 0, Am, = GeV, diag. 2-loop 
9 f = 0, Am, = GeV, RG 
9 r = -7t/4, Am, = 340 GeV, diag. 2-loop 
9 f = -it/4, Am, = 340 GeV, RG 



300 



500 



700 900 
m f [GeV] 



1100 1300 1500 



Fig. 4. Comparison between the Feynman-diagrammatic calculations and the re- 
sults obtained by renormalization group methods [8]. In the upper plot the (un- 
physical) soft SUSY-breaking parameters of the i-mixing matrix are chosen as 
input, while in the lower plot the physical i-masses and the mixing angle 9~ t are the 
input parameters. For the curves with #j = in the lower plot a mass difference 
Am,- = GeV is taken, whereas for 8^ = — 7r/4 we choose Am, = 340 GeV, for 
which the maximal Higgs-boson masses are achieved. 



5. Compact approximation formula for mh 



In order to extract the dominant contributions to from the rather 
complicated full result, we have derived by means of a Taylor expansion 
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a short analytical approximation formula from the diagrammatic two-loop 
result [16]. It can easily be implemented into existing programs and allows 
a very fast numerical evaluation. Since the most important contributions 
have been isolated in this analytical formula, it is also helpful for a better 
qualitative understanding of the source of the dominant corrections. 

In deriving the formula the following approximations have been made: 

• The momentum dependence of the one-loop and two-loop self-energies 
S s , s = h,H, hH has been neglected in eq. (1). 

• The parameters M, m g have been chosen according to M = nig = 
\J M| — to 2 , where Ms is given by 

M s = 



m| + mi : M h = M~ tR = m g ~ 

MlMf +mUMf + Mf) + mtY: M h * M iR 



Contributions from the t, t-sector up to the two-loop level: 

The main step of our approximations consists of a Taylor expansion 

of the one- loop and two-loop contributions from the t, t-sector in the 

parameter 

\m t M t LR \ = m| ~ ml 
Ml m~ + m 2 



A- - 1 — f - - Si (d) 

S> t2 tl 



where terms proportional to M| have been neglected in the t mass 
matrix. For the one-loop correction we have expanded up to 0(A~). 
We have kept terms up to 0{M^/mf), while terms of 0(M^/M^) 
have been neglected. For the two-loop self-energies the expansion has 
been carried out up to 0(A~). We have furthermore used the ap- 
proximation fj, = in the S s (0). After extracting a common prefactor 
(1/ sin 2 f3) we have set otherwise sin f3 = 1 in the non- logarithmic one- 
loop contributions, while the full dependence on sin [3 is kept in the 
logarithmic one-loop and the two- loop contributions. For a discussion 
of these approximations see Ref. [16]. 

For the one-loop contributions from the other sectors of the MSSM 
the leading logarithmic approximation has been used [5]. 

Corrections beyond 0(aa s ): 

Leading contributions beyond (D(aa s ) have been taken into account 
by incorporating the leading two-loop Yukawa correction of 0(Gj l m^) 
[7, 20] and by expressing the t, t-contributions through the MS top- 
quark mass m t instead of the pole mass m t according to eq. (3) . 
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The approximation formula for m\ is obtained by inserting the described 
approximations for the one-loop and two-loop self-energies S s into the mass 
matrix eq. (1). The diagonalization of the mass matrix incorporates con- 
tributions to m\ that are formally of higher order but are non-negligible in 
general. For large Ma these higher-order contributions are suppressed by 
inverse powers of Ma- Therefore it is possible for Ma S> Mz to perform an 
expansion in the loop order, leading to a very compact formula for m\ of 
the form 



m l = m 2 ' trec + ^.a.t/f + Am 2,a,rcst + Am 2,» s + ^ 



(5) 



The tree-level prediction and the one-loop contribution from the t, t-sector 
are given by 



m 



2,tree _ 



m| + m| 
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where A 



ifi 2 +^S 4 



3°w 



9°w . 



The dominant two- loop contribution of O(oa s ) to m\ reads: 



Am 



2,aa s 



GnV2a s a 
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,M t Lfi 
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For the one-loop contribution from the other sectors of the MSSM, Amjj' a ' rest . 

2 a 2 

and the leading two-loop Yukawa correction, Am h ' , which are numerically- 
less important than the contributions given above, we refer to Ref. [16]. In 
the contributions from the t, t-sector at one-loop and two-loop order, eqs. (7) 
and (8), we have included correction factors of 0{M^/M\). In this way the 
compact formula (5) gives a reliable approximation for Ma values down to 
at least Ma = 200 GeV. The approximation formula for the general case of 
Ma is given in Ref. [16] . 

The contribution of 0(aa s ) given in eq. (8) can be compared with ana- 
lytical formulas derived via the two-loop effective potential approach for the 
case of no mixing in the t sector [10] and via RG methods [7,9]. The leading 
term ~ log 2 (to 2 /M|) agrees with the results in Refs. [7,9, 10]. The sublead- 
ing term for vanishing t-mixing ~ log(m 2 /M|) agrees with the result of the 
two-loop effective potential approach [10] and the result of the two-loop RG 
calculation [9,10], but differs from the RG improved one-loop result [7,9]. 
The term ~ log(m 2 /Mj )(M/ y - R /Ms) 2 for non-vanishing i-mixing differs from 
the result given in Ref. [7,9]. All other terms of 0(aa s ) are new. The term 
~ M/^/Mg shows that the result for is not symmetric in ±Mj . The 
good numerical agreement with the RG results in the case of no mixing in 
the t sector can qualitatively be understood by noting that in the no-mixing 
case the leading term in both approaches agrees, while for the corrections 
proportional to powers of M^ R /Ms deviations occur already in the leading 
contribution. 

The compact approximation formula has been implemented into Feyn- 
Higgs, in order to allow a direct comparison between the full result and the 
approximation, and is also provided as a separate program called FeynHiggs- 
Fast. The improvement in the speed of the evaluation with FeynHiggsFast 
compared to FeynHiggs is about a factor of 3 x 10 4 . In Fig. 5 the approxima- 
tion formula is compared with the full result (in the lower plot the formula 
for general Ma is applied). The compact formula approximates the full re- 
sult better than about 2 GeV for most parts of the MSSM parameter space. 
Larger deviations can occur for \M^ R /rriq\ > 2. 

G.W. thanks M. Jezabek and the other organizers of the Epiphany Con- 
ference for the invitation, the excellent organization and the pleasant atmos- 
phere during the Conference. 
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